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I. INTRODUCTION 

The task of building a quantum computer, allowing 
for massive quantum parallelism pQ, long-lived quantum 
state superposition and entanglement, is one of the big 
challenge of twenty-first century physics. It requires pre- 
cise maintenance and manipulation of quantum systems 
at a level far from anything that has ever been done be- 
fore. Many criteria, nicely summarized by the DiVin- 
cenzo criteria [2 , need to be satisfied in order to be 
able to harness the power of quantum physics to ex- 
tract powerful information processing, and while most 
of these criteria have been satisfied to a reasonable level 
by some prospect for the implementation of a quantum 
information processing device [3 , no such prospect has 
yet been able to fulfill to a sufficient level all of these cri- 
teria at once. In particular, one criterion that not many 
prospects seem to be able to fulfill is that of scalability, 
the ability of the model to remain coherent for quantum 
systems with a large number of qubit s. 

Solid-state quantum registers seem to emerge as a 
medium that would allow for such scalability, while still 
being able to satisfy the other criteria [3]. However, as 
for any physical implementation of quantum computing, 
it suffers from noise, dissipation processes mainly in the 
form of decoherence and relaxation [4]. To be able to 
use the power of quantum computation, we must be able 
to eliminate this noise and maintain quantum informa- 
tion throughout the processing. A general framework for 
doing so is quantum error correction (QEC) jSJ [6]. It 
is known that by encoding a single logical qubit on suf- 
ficiently many physical qubits, we can correct arbitrary 
errors on a single physical qubit [5]. But what happens 
with more realistic multi-qubit errors? 

In order to perform a quantum error correction, ad- 
ditional qubits are needed. However, these additional 
qubits can lead to a much higher rate of decoherence 
(even superdecoherence) [Tj |8j . The implementation of 
the three-qubit repetition code has already been simu- 
lated for a cavity- QED setup [9], showing a significant 
but limited improvement for the preservation of quantum 



information. Here, we consider a realistic noise model for 
N charge qubit quantum registers [8], and simulate this 
noise model on encoded qubits to evaluate the efficiency 
of a QEC scheme. The quantum error correcting code 
(QECC) used is the 5 qubit DiVincenzo-Shor code [ID] . 

We structure this work as follows: we first describe the 
representation for quantum systems we will be using, as 
well as that for the noise model. We then further analyze 
the physical model we study, and review the QECC we 
implement within this model. Section |VI| then contains 
the core of this work, that is, a description of the algo- 
rithm used for the simulation, followed by a discussion of 
the results obtained. 



II. REPRESENTATION OF QUANTUM 
SYSTEMS 

We will be considering the basic unit of quantum infor- 
mation: the qubit. A qubit can be implemented by any 
two level quantum system, and so can be represented by a 
unit length vector in a two-dimensional complex Hilbert 
space, up to an equivalence relation for global phase. If 
we consider some prefered orthonormal basis states |0) 
and |1), a single qubit can be parametrized by 2 param- 
eters e [0,7r] and (j> e [0,2n): 

|^)=cosQ)|0)+e^sinQ)|l). (II.l) 

This is the Bloch sphere representation, for pure qubit 
systems. 

When considering multiple interacting qubit systems, 
the overall system is described by the tensor product of 
the component Hilbert spaces. In general, once two sys- 
tems have interacted, a complete description of the whole 
system cannot be obtained by a description of each sub- 
system: this is entanglement. Much of the power of quan- 
tum information processing comes from entanglement, it 
is a uniquely quantum resource. Note that this is also 
related to the difficulty of simulating a multi qubit quan- 
tum system: as a N qubit quantum system cannot be 
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completely described by each of the N two dimensional 
subsystems, it must be described by a 2^ dimensional 
system, so that the memory required to store the infor- 
mation about a N qubit quantum system scales exponen- 
tially with TV, and so do the time required to simulate 
operations on this system. 

This also raises the question of the description of sub- 
systems of entangled systems. It turns out that the den- 
sity operator formalism [TT] can solve this issue. In this 
formalism, a state ip is represented by a density matrix 
pi/j = I^X^I instead of by a state vector \ip). Then, if 
we are only interested in subsequent evolution of a par- 
ticular subsystem of the state -0, we can do a partial 
trace [11] over the other irrelevant subsytems and only 
evolve the reduced density matrix of the subsystem of 
interest. In this way, we obtain the right outcome statis- 
tic on this subsystem while not having to evolve a high- 
dimensionality state as required to describe the whole 
system completely. This density operator formalism also 
comes in handy for the description of statistical ensem- 
bles {pi,\^i)}, where the system is prepared in the state 
\ipi) with probability p^. The corresponding density op- 
erator is p — ^Pil^iX^il, and here also the evolution of 
this state leads to the right measurement statistics. 

A closely related subject is that of quantum noise. As 
a quantum system cannot ever be completely isolated 
from its environment, they interact and then become en- 
tangled. So, even though the global system-environment 
pair undergoes unitary evolution, the main system does 
not, and this generates quantum noise. This is well de- 
scribed by the quantum operation formalism [TT]. Two 
types of noise in which we will be mostly interested will 
be decoherence and relaxation, and we will see that these 
are closely related to generalized dephasing channels [12] 
and amplitude damping channels respectively [HI [13] . 

Since the quantum system of interest will be subject 
to noise, we want a way to quantitatively describe how 
close it will remain to the initial quantum information 
we wished to preserve. A good measure of this is the 
fidelity JTTJ [14], and for a pure input state evolving 
to a mixed state p, the fidelity is given by 



F(M,p) = vW). (n.2) 



III. QUANTUM NOISE 

The two main types of noise in which we will be in- 
terested are decoherence, which correspond to a loss of 
quantum coherence without loss of energy, and relax- 
ation, which correspond to a loss of energy. We will 
see that these types of noise are closely related to well- 
studied channels: generalized dephasing channels [T2| fT5] 
and amplitude damping channels JTTJ [13] 



A. Generalized dephasing channels 

Generalized dephasing channels correspond to physi- 
cal processes in which there is loss of quantum coherence 
without loss of energy. That is, there exist a preferred 
basis, called the dephasing basis, such that pure states 
in that basis are transmitted without error, but pure su- 
perposition get mixed, and quantum information is lost 
to the environment. 

If we let A be the input system with orthonormal basis 
B be the output system with orthonormal basis 
{\i) B }, and E be the environment with normalized, not 
necessarily orthogonal states {|<^) }, then an isometric 
map from the input system to the output-environment 
system is given by 

u A ^ BE = ^2\i) B m E (i\ A . 

Then, starting with an input state p A , we get the output 
state a B by first applying the isometry JJ A ^ BE ', then 
tracing over E: 

a B = Tr E (Up A tf) = ^> | / \J) l*)01 B &l&>- (IH-2) 

As we can see, in this representation, the output corre- 
sponds to the Hadamard product between the ((i | p A \j)) 
input matrix and the decoherence matrix. Also, 

since the |^) are normalized, = 1, the diagonal 

terms are left unchanged: 

(i\a B \i) = (i\p A \i). (III.3) 

These diagonal terms correspond to the probability of 
measuring the quantum state in the corresponding de- 
phasing basis state, while the off-diagonal terms are co- 
herence (phase) terms. For these off-diagonal terms, the 
Cauchy-Schwarz inequality (|(&|£j)| 2 < = 
1) tells us that the output terms must be smaller than 
or equal to the input terms: | (i | a B \j) | < | (i | p A \j) |. 
In the special case where the {1^)^} are also orthogonal, 
we obtain a completely dephasing channel, which kills 
all off-diagonal terms, (i | a B \j) = if i ^ j. For a sub- 
sequent measurement in the dephasing basis, this yields 
a classical statistical output with no quantum coherence 
terms: a B I p A \i) \i)(i\ B . 

B. Amplitude damping channels 

Amplitude damping channels correspond to physical 
processes in which there is a loss of energy from the 
quantum system of interest to the environment. That is, 
there is a probability that a quantum state passes from a 
higher energy excited state to a lower energy state. For 
the two-dimensional quantum systems in which we are 
mostly interested, this corresponds to a probability of 



passing from the excited state to the ground state, while 
the ground state is left unaffected. In the quantum op- 
eration formalism, this channel has operation elements 
{£o,£i}: 



-En = 
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(III.4) 



where the matrix representation is given in the energy 
eigenbasis, and the parameter 7 corresponds to the prob- 
ability of decay. Then, for some input state p A , we get 
an output state a B : 



£ E ^ AE l 



(III.5) 



=0,1 



If we represent the action of the channel as £ , i.e. 
£(p A ) = cr B , we can extend this definition to multi-qubit 
systems. Making the assumption that energy lost for 
each two-dimensional subsystem is independent, a ten- 
sor product of this channel, £® n , corresponds to a loss of 
energy in each subsystem independently. 




FIG. 1: Physical model for the 5-qubit quantum register. The 
qubit geometry is inspired from recent experiments on coher- 
ent level mixing in vertical quantum dots [22] . 



This leads to the decay of the off-diagonal elements of 
the N-qubit density matrix. The structure of this decay 
is non-trivial and is shown in FIG. [2] In general, this 
leads to an error, which cannot be reduced to a single 
qubit error. 



IV. DECOHERENCE AND RELAXATION IN 
CHARGE QUBIT REGISTERS 



Solid state quantum systems represent an interest- 
ing prospect for the physical implementation of scalable 
quantum devices [3 . A lot of research goes into imple- 
mentation of superconducting qubits [161 [17] . Another 
important avenue of research are quantum dots, in which 
two level quantum systems can be implemented in either 
the spin or position degree of freedom of the electron 
trapped in the dots. While spin quantum dots have a 
longer decoherence time [18j [19] , spins are also more dif- 
ficult to manipulate. In charge quantum dots, where the 
basis states correspond to the position of the electron 
in either of two adjacent quantum dots, the coupling is 
strong and hence fast manipulation is possible, but this 
also leads to shorter decoherence times. We will consider 
such a system, which offers an interesting playground for 
current technology research [20j[2T]. Moreover, the de- 
coherence for N qubit quantum registers, which can be 
seen in FIG. [I] has already been computed by Ischii et al. 
[8], and it is this model which is used in this simulation. 

It is found that for these charge quantum dots, the de- 
coherence of the quantum system can be represented as 
a generalized dephasing channel in the canonical basis of 
the charge qubits, corresponding to the position eigen- 
states. Indeed, pure states in that basis are perfectly 
transmitted, while pure superpositions become noisy. It 
is then possible to compute the evolution of the quantum 
register by taking the Hadamard product of the input 
density matrix with a decoherence matrix computed in 
[8], thus getting the decohered output state. We chose 
as physical timescale unit c^o, which corresponds to the 
upper cut-off frequency. In charge qubits in GaAs we 
typically have 

uj = 0.2 x 10" 10 8. (IV.l) 
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FIG. 2: (Color online) Matrix elements of the decay function 
density matrix due to decoherence. 

In addition, this decoherence model does not take into 
account any relaxation the system might be subject to. 
We assume that the relaxation is independent of deco- 
herence, and also that energy loss for each qubit is inde- 
pendent. For each two dimensional quantum subsystem, 
if |0) and |1) correspond to the two position eigenstates, 
then the two energy eigenstates are |+) = 77^ (|0) + |1)) 

for the ground state, and | — ) = -^=(|0) — |1)) for the 
excited state. We can then compute the effect of relax- 
ation by operating a tensor product of amplitude damp- 
ing channels to e ach qu bit in the register, where the op- 
eration elements (III. 4) are in the {|+), |— )} basis. The 
decay parameter 7 is time (t) and temperature (T) de- 
pendent, 



7 = 1 — exp(—t • T), 
setting ks = h = 1. 



(IV.2) 
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V. QUANTUM ERROR CORRECTION WITH A 
PERFECT 5 QUBIT CODE 

Since our quantum register is subject to noise, we want 
to encode an input quantum state in such a way that 
we can preserve the quantum information even when the 
physical system undergoes noise. A general framework 
for doing this is quantum error correction (QEC) [5j|6], 
which encode a logical qubit containing the quantum in- 
formation on multiple physical qubits, and then usually 
by performing syndrome measurements on these physi- 
cal qubits, we can determine a certain set of errors, and 
apply the corresponding correction [IT]. Many quantum 
error correcting codes (QECC) have been found that can 
correct an arbitrary error on a single physical qubit, and 
it is known that the smallest such codes require 5 physical 
qubits [23]. 

One of these 5 qubit codes is the DiVincenzo-Shor 
QECC [10]. The code words for the logical states 
{|0) L ,|l) L }are [3: 

|0) L = -( 100000) + |11000) + |01100) + |00110) 

+ |00011) + |10001) - |10100) - |01010) 

- |00101) - |10010) - |01001) - |11110) 

- |01111) - |10111) - |11011) - |11101) ) 

|i) L = i( |imi) + loom) + |iooii) + |iiooi) 

+ |11100) + |01110) - |01011) - |10101) 

- |11010)-|01101)-|10110) 

- |10000) — |01000) — |00100) ■ 
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with stabilizers 



M = I®Z®X®X®Z 
M 1 = Z®I®Z®X®X 
M 2 =X(g)Z(g)/(g)Z(g)X 
M3 = X®X®Z®I®Z, 



(V.l) 



which leave the code invariant. The X and Z opera- 
tors are respectively the X and Z Pauli operators. The 
encoding circuit can be seen in [3 . This is a perfect 
non-degenerate code, meaning that all single qubit er- 
rors map to a different syndrome, as can be seen in TA- 
BLE [I] To obtain a measurement of the error syndrome, 
we measure observables in the eigenbasis of each of the 
four stabilizers in ( |V.1[ ), and depending on the outcome 
of these four measurements we apply the corresponding 
correction from TABLE [H 

An alternative to these multi-qubit measurements [11 
is to use ancillary qubits to conditionnally apply each 
of the four stabilizers on a |+) state, thus recording the 
phase (M0-M3 have eigenvalues ±1). Measuring these 
four ancilla qubits in the {|+), | — )} basis then gives the 
same result as above without having to perform multi- 
qubit measurements. This syndrome measurement cir- 
cuit is given in [3]. 



TABLE I: Table of the different possible single qubit errors 
Em, with their associated syndrome measurement output M, 
which has binary expension M3M2M1M0. 



This code is guaranteed to correct any single qubit 
error, but it is not designed to correct multi-qubit errors, 
which our noise model will introduce. It will then be 
interesting to see how well this code will perform for an 
approximate correction of multi-qubit errors. 



VI. QUANTUM ERROR CORRECTION IN A 
CHARGE QUBIT QUANTUM COMPUTER 

The qubit model for physical implementation we con- 
sider in this simulation is a charge quantum dot. The two 
principal sources of noise with this model are decoherence 
and relaxation. 

The effect of decoherence on a iV-qubit register had 
been computed in [8], where decoherence is caracteristic 
of a generalized dephasing channel [12], leaving compu- 
tational basis states unaffected. 

The effect of relaxation correspond to a probabilistic 
loss of energy, where the excited state correspond to the 
|— ) state, and the ground state correspond to the |+) 
state. We make the assumption that this relaxation for 
each qubit is independant, so that we implement this 
noise has a tensor product of an amplitude damping 
channel in the {|+), |— )} basis, on each qubit. 

To counter the effect of noise, we encode a logical qubit 
into 5 qubits, using the DiVincenzo-Shor QECC, which 
can perfectly correct errors on a single qubit. However, 
the noise model used creates errors on multiple qubits, 
and so it is not obvious that QEC with the DiVincenzo- 
Shor code can extend the lifetime of our logical qubit. 

Here, we verify the effect of such a quantum error cor- 
rection scheme, for different input states and different 
frequencies of correction. The algorithm used for the 
simulation is presented first, and the results of the simu- 
lation are presented next. 



A. Implementation of QEC algorithm 

An overview of the algorithm can be seen in FIG. [3] 
First, the logical qubit is encoded into 5 physical qubits 
using the DiVincenzo-Shor QECC. This encoding is as- 
sumed to be perfect. We only consider pure input qubits, 
which are completely caracterized by two real parame- 
ters corresponding to their position on the Bloch sphere: 
e [0,tt], </) e [0,2tt). 

After encoding, the quantum register is subjected to 
noise for a given period of time t, depending on the fre- 
quency of error correction. In our simulation, the de- 
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where the order in the operator product is irrelevant since 
the terms commute. 

Following the quantum error correction step, a record 
of the fidelity of this output state compared to the input 
state is taken. We then reloop over the noise, quantum 
error correction and fidelity recording steps, until the sum 
of all the noise time steps add up to the desired total time 
of the simulation. 



FIG. 3: Overview of the algorithm for the simulation of quan- 
tum error correction of a solid-state quantum register. 



coherence noise and the relaxation noise are applied suc- 
cessively, and we checked numerically that the order does 
not affect the result. The intensity of the noise depends 
on the temperature T of simulation, and the probability 
of relaxation is adjusted such that without error correc- 
tion, after a time ujo the fidelity of the noisy state with 
respect to the pure input state is the same as that of 
decoherence. 

Quantum error correction is then performed. Here we 
assume perfect operation of each component of the QEC 
algorithm. However, instead of performing the measure- 
ment of the syndrome and then applying the correspond- 
ing correction (see TABLE [i]), we defer the measurement 
to the end of the circuit and perform conditional quan- 
tum correction [3l [TT] . To do so, four ancillary qubits are 
used to record the syndrome measurements. It is pos- 
sible to avoid multi-qubit measurements by using these 
ancillary qubits, but the drawback is that we pass from 
a 5 qubit quantum system to a 9 qubit quantum system. 
This increases the difficulty of physical implementation 
because of the additionnal qubits, and also because every 
round of QEC requires four fresh ancilla qubits. This also 
slows the simulation a lot, since we pass from a 2 5 = 32 
dimensional complex vector space to a 2 9 = 512 one, and 
the matrix operations are accordingly scaled. 

After performing the conditional QEC, these ancilla 
qubits are traced over, which corresponds to a measure- 
ment without a record of the outcome. This gives an 
output state which is a statistical ensemble correspond- 
ing to the different possible corrected states. 

To perform the conditional quantum error correction, 
we implement a quantum circuit which would correct for 
the different errors in TABLE [TJ corresponding to the 
value of the ancilla qubits. If we let \M) correspond to 
the ancilla state with binary expansion M3M2M1M0, for 
M running from to 15, and if Em is the corresponding 



B. QEC performance 

1. Decoherence 

We can see on FIG. 0] the results of the simulation 
when our quantum register is subject only to decoher- 
ence. With a small frequency of error correction of 
co^ 1 , the fidelity is already improved over the uncorrected 
qubit. For a basis qubit, either |0) or |1), we do not have 
decoherence, but the encoded state does suffer from de- 
coherence. The encoded basis states however do exhibit 
an interesting behavior when they are decohered, but not 
error corrected. In fact, even though the noise takes the 
5 qubit state out of the code, when decoded we still get 
back the original basis state. For these basis states, we 
could understand this behavior by noting that the code 
words for the logical |0) and | 1) con tain different 5 qubit 
basis states, as we can see in (V.l). But for pure super- 



position of these basis states, if we do not error correct 
the encoded states, it also decoheres in a way such that 
once decoded, the same state as the unencoded decohered 
state is obtained, which is even more surprising. 



2. Relaxation 

We can see in FIG.[5]the results of the simulation when 
we only consider relaxation affecting our quantum regis- 
ter. This type of noise has a stronger temperature depen- 
dency and is stronger than decoherence for temperatures 
above 5 mK, at which point we obtain on average a sim- 
ilar drop in fidelity at 1 ujq. It is also harder to correct 
than decoherence. We thus need a much higher frequency 
of quantum error correction to keep the fidelity close to 
1 until we have reached a saturation in relaxation, that 
is we have reached the lowest energy state for the un- 
encoded qubit. This low energy state is the only state 
invariant under action of the channel. 
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Decoherence only (T=5mK) 




time [©J 

FIG. 4: Simulation of quantum error correction with only 
decoherence. The time between two successive QECs is 
At = cjo- For a total of IOOOcjo, the top curve shows the 
result of the QEC and the non-corrected lqubit decoher- 
ence is shown in the bottom curve. The input qubit is 
= cos (|) |0)+e*sin(!) |1) and the temperature is as- 
sumed to be 5 mK. 



Relaxation only (T=5mK) Relaxation only (T=5mK) 
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FIG. 5: Simulation of quantum error correction with only 
relaxation: In the left graph the fidelity as a function of time 
is shown for the same At = cjo as in FIG. [4] By reducing the 
time interval between two QECs the fidelity remains close to 
one as seen in the right graph. For both graphs, the input 
qubit is the same as in FIG.^] 

3. Simultaneous decoherence and relaxation 

When combining the effect of decoherence and relax- 
ation, the frequency of error correction is on the order of 
what is needed for relaxation, since relaxation is much 
harder to correct than decoherence. Here too, when tak- 
ing into account decoherence and relaxation, we are still 
able to maintain the fidelity arbitrarily close to one (see 
FIG.[6|. This tells us that even though it may require a 
high frequency of quantum error correction, the 5 qubit 
code is able to correct a realistic multi-qubit error sys- 
tem. 



C. Discussion 

We have only considered noise in the N-qubit system. 
That is, an important assumption made for our simula- 
tion is that of a perfect operation of the QEC algorithm. 




time [©J 

FIG. 6: 1 minus fidelity is shown as a function of time for a 
very short At = O.OOIcjo- The high-frequency QEC for joint 
decoherence and relaxation enables the quantum information 
to be preserved almost perfectly. The input qubit is the same 
as in FIG. H 

However, in a realistic fault-tolerant quantum compu- 
tation (FTQC) [25] , errors will also occur within the 
operation of the gates and measurements, and this also 
during the QEC steps. Similar concerns led some physi- 
cists to wonder if the error model adopted to get to the 
threshold theorem [26] [27] for FTQC is physical enough 
so that FTQC is possible at all [28]. Hence, such a simu- 
lation, with a realistic model for operations, would be an 
important step toward determining if physical implemen- 
tations of quantum computers are possible. It is possible 
that these types of errors can be corrected in a similar 
manner, but this is beyond the scope of this work. 

The model considered introduces multi-qubit errors, 
which are much smaller than full single-qubit ones when 
the time evolution is very short. At the lowest order, 
we showed that these small multi-qubit errors can be 
corrected by single-qubit QEC codes. However, this in- 
creases the required QEC frequency, but it is certainly 
possible to find other error correcting schemes optimized 
for multi-qubit errors, which would require a lower QEC 
frequency. 

Also, since the decoherence time for charge qubit s is 
short, the frequency of QEC required to maintain good 
fidelity for extended periods is too high to be feasible 
experimentally with current techniques. However, this is 
a proof of principle and it is reasonable to expect that 
for other implementations of qubits like spin qubits or 
superconducting qubits, where decoherence times can be 
much longer, similar performances would require a QEC 
frequency that is attainable experimentally. 



VII. CONCLUSION 

We have simulated the implementation of the 5 qubit 
DiVincenzo-Shor QECC for a charge qubit quantum reg- 
ister. The register was subjected to a realistic noise 
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model, consisting of decoherence and relaxation. Even 
though the QECC is designed to perfectly correct only 
single qubit errors, an application of the QEC routine 
with a high enough frequency was able to limit the effect 
of multi-qubit errors introduced by our noise model. In 
fact, by adjusting the frequency of error correction, we 
can get the fidelity as close to one as possible, and for 
extended periods of time. When considering both types 
of noise separately, relaxation showed to be harder to er- 
ror correct than decoherence, requiring higher frequency 



of QEC for the same results. 
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